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ABSTRACT 

Based on the diffusion-halo model for cosmic-ray (CR) propagation, includ- 
ing stochastic reacceleration due to collisions with hydromagnetic turbulence, 
we study the behavior of the electron component and the diffuse 7-rays (D7's) 
induced by them. The galactic parameters appearing in these studies are es- 
sentially the same as those appearing in the hadronic CR components, while 
we additionally need information on the interstellar radiation field, taking into 
account dependences on both the photon energy, E p ^ and the position, r. We 
compare our numerical results with the data on hadrons, electrons and D7's, 
including the most recent results from FERMI, which gives two remarkable re- 
sults; 1) the electron spectrum falls with energy as E~ 3 up to 1 TeV, and does not 
exhibit prominent spectral features around 500 GeV, in contrast to the dramatic 
excess appearing in both ATIC and PPB-BETS spectra, and 2) the EGRET 
GeV-excess in the D7 spectrum is due neither to an astronomical origin (much 
harder CR spectrum in the galactic center) nor a cosmological one (dark matter 
annihilation or decay), but due to an instrumental problem. In the present paper, 
however, we focus our interest rather conservatively upon the internal relation 
between these three components, using common galactic parameters. We find 
that they are in reasonable harmony with each other within both the theoretical 
and experimental uncertainties, apart from the electron-anomaly problem, while 
some enhancement of D7's appears in the high galactic latitude with |6| > 60° in 
the GeV region. 

Subject headings: cosmic rays — Galaxy: structure — electrons: diffuse back- 
ground 
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1. Introduction 

Although the electron component is only a small fraction of all cosmic-ray (CR) compo- 
nents, around 1% of the proton intensity around 10 GeV, it plays a key role in understanding 
the structure of our Galaxy and the galactic phenomena occurring within it. This is because 
electrons have electromagnetic interactions with the interstellar radiation field, such as pho- 
tons and magnetic fields, resulting in drastic energy loss during propagation through the 
Galaxy, in contrast to the hadronic component. 

This peculiar nature yields valuable information for the study of CR astrophysics, which 
can not be obtained by the hadronic components alone. Namely, due to the rapid energy-loss 
rate, proportional to E\ in the high energy region, from the inverse Compton scattering off 
photons and synchrotron radiation in magnetic fields, the life-time of TeV electrons is at most 
10 5 yr, indicating that detected electrons have originated in nearby sources, less than 1 kpc 
from the solar system (SS). Therefore, accurate observations of TeV electrons will provide 
a direct signature of nearby CR sources as well as the mechanism of the CR acceleration, 
while depending on the release time from supernova remnants and their distance from the 
SS. 

Qualitative studies of such a possibility have been performed by many authors (Shen 
1970; Nishimura et al. 1979; Cowsik & Lee 1979; Berezinskii et al. 1990; Aharonian et al. 
1995; Ptuskin & Ormes 1995; Pohl k Esposito 1998; Kobayashi et al. 2004; Delahaye et 
al. 2010), with Kobayashi et al. and Delahaye et al. presenting explicitly several candidates 
for nearby sources of high energy CR electrons, based on the most recent data for the age 
and distance of each supernova remnant near the SS, although the statistics of high energy 
electron data are currently too poor to identify sources definitely. 

Particle identification and the energy determination of high energy electrons is, however, 
quite difficult, while direct observation of low energy electrons is relatively easy using, for 
instance, magnetic spectrometers, and has been performed by several groups (Golden et al. 
1994; Boezio et al. 2000; DuVernois et al. 2001; Aguilar et al. 2002). 

Although the statistics are not sufficient, the only group that succeeded in observing 
directly TeV electrons is Nishimura et al. (1980; see also Kobayashi et al. 1999) with the use 
of the balloon-borne emulsion chamber. It should be noted that they actually observe event 
by event the vertex point of the electron with subsequent e ± -pair due to bremsstrahlung 7, 
with no uncertainty from proton contamination. The precision in the energy determination 
is approximately 10% for electrons in the energy region larger than 50 GeV, based on both 
the three-dimensional cascade theory (Nishimura 1964) and the simulations (Kasahara 1985; 
Okamoto & Shibata 1987), which have been well established by the use of accelerator beams 
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(Hotta et al. 1980; Sato & Sugimoto 1979). 

Recent development in high energy electron observations is indeed remarkable, partic- 
ularly those of ATIC (Chang et al. 2008) and PPB-BETS (Torii et al. 2006), which showed 
an anomaly in the electron spectrum with a significant bump around 500 GeV. Both groups 
point out that the excess indicates either a nearby source of energetic electrons, or those 
coming from the annihilation of dark matter particles. 

On the other hand, the most recent results obtained by the FERMI Large Area Telescope 
(FERMI-LAT; Abdo et al. 2009) present no prominent excess, with the electron spectrum 
falling with energy as E~ 3M up to 1 TeV, which is not inconsistent with the emulsion cham- 
ber data (Kobayashi et al. 1999) within the statistical errors. The H.E.S.S. ground-based 
telescope (Aharonian et al. 2008, 2009) also shows no indication of structure in the electron 
spectrum, but rather a power-law spectrum with ^- 30±0 - 1±0 - 3 (o.l: stat. error, 0.3: syst. 
error), albeit this being an indirect observation. 

Nevertheless, looking carefully FERMI data around the anomaly-energy, they still show 
systematically an enhancement as large as 30% compared to the numerical results (Abdo et 
al. 2009; Strong et al. 2004; see also Figure 14 in this paper), so that we can not exclude 
the possibility of an additional component such as local sources and/or the dark matter 
scenario, while strength of the anomaly compared to the background diffuse electrons is not 
as dramatic as presented by ATIC and PPB-BETS. 

In any case, both observational and theoretical studies for high energy electrons are 
becoming increasingly important not only for astrophysics, but also for particle physics and 
cosmology. It is, therefore, desirable to find a reasonable model for electron propagation 
in the Galaxy, which must explain consistently and simultaneously all CR observables and 
not just electrons, using common galactic parameters with the smallest number of variables 
possible. In the sense, the recent review article by Strong, Moskalenko, & Ptuskin (2007) 
is a useful survey of both the theory and relevant experimental data for the propagation of 
CRs, comprehensively summarizing the current landscape and open questions, although it 
was published just before the anomaly problem mentioned above. 

Under these situations, we have studied the three-dimensional CR propagation model 
analytically, and found excellent agreement with the experimental data for various hadronic 
components, stable primaries, secondaries such as boron and sub-iron elements (Z =21-23), 
isotopes such as 10 Be, and antiprotons as well, in four papers, (Shibata et al. 2004, 2006, 
2007a, 2008), hereafter referred to as Papers I, II, III and IV, respectively. 

We have applied our model further to the studies of diffuse 7-rays (D7's) (Shibata, 
Honda, & Watanabe 2007b; hereafter Paper V), and found that all these components are 
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generally in agreement with each other using the same galactic parameters, within the uncer- 
tainties in the experimental data and various kinds of cross-sections used for the numerical 
calculations. However, in Paper V, we use the simulation results for electron-induced 7-rays 
provided by Hunter et al. (Bertsch et al. 1993; Hunter et al. 1997), where the modeling 
of CR propagation and the galactic parameters assumed are somewhat different from ours. 
So we have yet to see complete internal consistency among all CR components — hadrons, 
electrons and D7's — using the same galactic parameters in our propagation model. 

In the present paper, we extend it to the electron component, based on the diffusion- halo 
model proposed by Ginzburg, Khazan & Ptuskin (1980), taking the reacceleration process 
into account. However, we focus in the present work on diffuse electrons in the steady 
state without discriminating those produced by nearby sources from those of distant ones, 
and present the intensity of the D7's produced by them in the energy range, E 1 = 30 MeV- 
100 GeV, covered by EGRET and FERMI. Comparison with radio and TeV-7 data will be 
reported separately in the near future. 

In order to apply our model to the electron component and electron-induced D7's, we 
need information on the interstellar radiation field (ISRF) in addition to the interstellar mat- 
ter (ISM), particularly their spatial gradients for the study of the (/, ^-distribution of D7's 
(I: galactic longitude; b: galactic latitude). Nowadays the most advanced and standard code 
for the ISM and ISRF models is GALPROP, extensively developed by Strong & Moskalenko 
(1998), incorporating the latest survey data in the very wide wavelength range from ultra- 
violet to radio. In the present work, we assume empirical density distributions for the ISM 
and ISRF, smoothing the numerical data given by GALPROP available most recently (Porter 
et al. 2008), in order to combine with our analytical solution for electron- induced D7's. 

In §2, we discuss the interstellar environment provided by GALPROP, focusing on 
the spatial distribution of both matter (atomic, molecular, and ionized hydrogen) and pho- 
tons (ultraviolet, visible, infrared, mid- and far-infrared, and cosmic microwave background 
[CMB] radiations), and in §3 we present the relevant elementary processes for electrons, fo- 
cussing on the energy losses due to ionization, bremsstrahlung, synchrotron, inverse Compton 
(IC), and on the energy gain due to the reacceleration. 

In §4, we present the diffusion equation, and give its solution explicitly in the steady 
state, N e (r;E e ), where the Klein-Nishina effect is quite important in the electron energy 
spectrum in the high energy region, > 10 GeV. In § 5, we present the emissivity of electron- 
induced 7's, g 7 (r; E~f), with use of realistic spatial distributions of ISM and ISRF as discussed 
in §2, and show the numerical results at several observational points, where those of the 
hadron-induced 7's are presented as well. In §6, we first summarize the galactic parameters 
and their explicit values expected from the CR data, and then compare our numerical results 
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of electron flux and D7's with recent observational data, including those most recently ob- 
tained by FERMI and H.E.S.S. Finally in § 7, we summarize the results, and discuss several 
remaining open questions, while we do not touch upon the so called electron-anomaly. 



2. Interstellar environment of our Galaxy 

2.1. Interstellar matter 

First we consider the ISM for two processes, ionization and bremsstrahlung. In Figure 1 
we plot histograms of column density for HI and H2 in the galactic plane (GP) given by 
GALPROP, where we also plot the empirical curves used in the present work, 
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with r in kpc, and p h ("/i" = HI, H 2 ) in 10 20 H atoms cm" 2 . The numerical values of the 
coefficients are summarized in Table 1. However, the choice of above empirical form is not 
critical, and other choices may be possible. 

The H 2 gas is strongly confined to the GP and its vertical structure is modeled by a 
gaussian distribution with a width of approximately 70 pc, while the HI gas lies in a flat 
layer with a FWHM of 230 pc in 3.5 kpc < r < r Q (=8.5 kpc), and is approximated by the 
sum of two gaussians and an exponential tail (Ferriere 2001; Moskalenko et al. 2002). Taking 
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Fig. 1.— Column density of interstellar hydrogen. Curves are empirical ones given by 
equation (1) with the parameterization summarized in Table 1. 
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Table 1. Summary of the numerical values of the coefficients appearing in equations (la) 
and (lb) in units of 10 20 H atoms cm~ 2 , where "(±m)" denotes the multiplication of 10 ±m . 
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Table 2. Summary of functions for ISM gas density, HI and H 2 , where r, r (=8.5kpc), z, 
and 2 are all in units of kpc, and nf in units of H atoms cm' 3 . 
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these situations into account, we assume the following spatial distribution for the ISM gas 
density, corresponding to equations (la) and (lb), 

2t£L ff 4 g£L , C*^HI,H 2 ), (2) 
< ~h(r©, 0) p h (r ) 

where n® (n® 2 ) is the gas density of HI (H 2 ) at the SS with typically n® w n® ps 0.5 
H atoms cm -3 . See Table 2 for the explicit forms of S m and E H2 . 

For the ionized hydrogen gas, nil, we use the two-component model of Cordes et al. 
(1991), n mi (r) = n^^r) +n^l(r), and both components are modeled by a gaussian-type 
distribution for the radial structure, and by a simple exponential one for the vertical struc- 
ture. The explicit values of the two components at the SS, [n^^o), n^r©)], are [0.025, 
0.013] cm -3 respectively (Cordes et al. 1991; Strong et al. 1998). So the contribution of HII 
is much smaller than those of H I and H2 and is not important in the present work. 



2.2. Interstellar radiation field 

First we consider the medium — virtual photons induced by the static magnetic field 
- for the synchrotron process. It is approximately given by an exponential- type gradient, 
while the scale height is not yet clear. Practically, for the study of synchrotron radiation, 
we need the energy density of virtual photons at r, e B (r), and assume in the present work 

e B (r*) = e B , exp[-(r/r B + \z\/z B )], (3) 

with 

e B ,o = B^/8-k, 

where B is the magnetic field at the galactic center (GC), and e B ,o is its energy density, for 
instance e B ,o ~ leVcm~ 3 for B = QfiG, and typically [2r B , 2z B ] [10, 2] kpc (Strong et al. 
2000). 

On the other hand, the photon gas for the IC process is somewhat different from those 
discussed above. Namely, we need the number density of the photon gas in the ISRF, 
n p h(r; .Eph), as a function of the target photon energy at r. Separating it into two 
parts, a r-dependent energy-density term, e p h(r*)> an d a r-independent term, W^ik) with 
k = E ph /[k B T ph ], we rewrite n ph (r; E ph ) as 
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E ph n ph (r; E ph )dE ph = e ph (r)W ph (k)d In k, (4) 

where k B is the Boltzmann constant, and T ph is the characteristic temperature of the ISRF. 

There are three main radiation sources in the photon gas, (i) the 2.7 K CMB radiation, 
(ii) stellar radiation with wavelengths of 0.1-10 /zm (ultraviolet-visible-near-infrared), and 
(hi) re-emitted radiation from dust grains at 10-1000 /im (mid-to-far-infrared). 

We classify them further into six wavelength bands, each labeled with i — for (i), % — 1, 
2, 3 for stellar-1, -2, -3 in (ii), and i = 4, 5 for dust-1, -2 in (iii) (see Fig. 2). Needless to 
say, there is no spatial gradient in the CMB (i = 0), which is distributed uniformly in space, 
e p h(r) = e p k = 0.261 eVcm -3 , and the normalized spectrum, W^(k), is given by the familiar 
Planck formula with T ph = 2.73 K. 

On the other hand, for (ii) and (iii) in the wavelength range A = 0.1-1000 /zm, the energy 
density, e p l(r) (i = 1-5), must depend on r, and wQ(k) is unlike the simple CMB spectrum, 
and is very complicated. In the following discussions, we often omit the suffix i for simplicity 
unless otherwise specified. 

In the present work, we assume a gaussian-type distribution in In A; for W p h(k), 

W ph (k) = -L- e -( lnfc ) 2 /^ 2 ); h = A /A, (5) 
V2na 



(6) 



so that the mean radiation intensity, J p h, is given by 

[m(A /A)P 



47rA/ ph (r; A) 1 

exp 



2a 2 



ce ph (r) V2tt<j 
where A is the peak wavelength for each radiation with k B T ph = 2nch/\ Q . 

In Figure 2, we present examples of the mean radiation intensity (multiplied by 47rA) for 
the maximal metalicity gradient {filled symbols) and no metalicity gradient (open symbols) 
at two radial distances, r = (squares) and 8kpc (circles) in the GP given by GALPROP, 
where also drawn are curves expected from the right-hand side of equation (6) for W p h(k), 
assuming 

e P h(?") = e P h,o exp [-(r/r ph + \z\/z ph )], (7) 

for r > 3kpc, and see caption of Table 3 otherwise. In Table 3 we summarize numerical 
values of [Xq \ T h af, e p ^ ] (2 = 0-5) with r ph = 3.2kpc irrespective of the population i, and 
also presented are those of e p l(r) for r < 3kpc, while they are independent of r except e p ^(r). 

Let us demonstrate the energy density separately for the stellar and the dust radiation, 
X)i=i e ph an d Si=4 e ph respectively, against the galactocentric distance r in Figure 3, where 



- 9 - 




wavelength; Z(jim) 

Fig. 2. — Interstellar radiation field (ISRF) at two galactocentric distances obtained by 
GALPROP, r = (GC; square symbols) and 8kpc (near SS; circle symbols). Open marks 
correspond to maximum metalicity gradient, and filled ones to the minimum metalicity 
gradient. Dotted curves are given by equation (6) with parameters summarized in Table 3 
for each population i, while the solid ones are those superposing them. CMB radiation (solid 
curve) is also shown for reference. 

we plot also numerical data given by Mathis et al. (1983; filled grey symbols). Two curves for 
the stellar emission and the dust re-emission are drawn by the use of the parameterization 
summarized in Table 3, where we do not take the difference in the choice of metalicity 
gradient into account, as it is effective only near the GC for the dust re-emission and is 
approximately one order of magnitude smaller than the stellar radiation. 

For the latitudinal scale height, z p ^, in equation (7), we assume z p ^ ~r p h/8 = 0.4kpc, 
referring to the speculation by Freudenreich (1998) based on the DIRBE (Diffused Infrared 
Background Experiment) survey, while the surveys of the diffuse FIR/sub-mm emission for 
the latitudinal direction at various radial distances r are not sufficient to construct a reliable 
model. 

3. Energy loss and gain 

3.1. Energy loss in ISM and ISRF 

The energy loss processes for the electron component are dramatically different from 
those for the hadronic components, with four main processes: bremsstrahlung (= "rad"), 
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Fig. 3. — ISRF energy density as a function of galactocentric distance r at z = for the 
stellar radiation and the re-emission from dust grains, each with the maximum metalicity 
gradient (filled black circles) and minimum metalicity gradient (open circles), where also 
shown are those obtained by Mathis et al. (1983) (filled grey circles). Solid curves are the 
empirical ones obtained by equation (7). 



ionization (= "ion"), synchrotron and IC (together = "sic"). For the bremsstrahlung (Koch 
& Motz 1959; Gould 1969; Ginzburg 1979), 



1 I AE P 



E e \ At , ' rad 



n(r)w Tetd (E e 



1+0 



n 



with 



n(r) = n m (r) + n HU (r) + n H2 (r), 



(9) 



where w ra d(E e S> m e c 2 ) = ~ 7.30 x 10 16 cm 3 s 1 , independent of E e with the complete 
screening cross-section in the high energy region; see Appendix A for the explicit forms of 
w vad (E e ), and §2.1 for n h , ("h" = HI, HII, H 2 ). 



Similarly for the ionization, we use the Bethe-Bloch formula (Ginzburg 1979) 
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Fig. 4. — Energy losses per unit time of CR electrons in ISRF at SS (r Q = 8.5kpc), shown 
separately for four components, synchrotron (heavy dotted line), IC's by stellar radiation 
(broken curve), by re-emission from dust grains (dotted curve), and by CMB (broken dotted 
curve), together with the sums, IC-all = IC-stellar + IC-dust + IC-CMB (thin solid curve), 
and synchrotron + IC-all (heavy solid curve). 



On the other hand, the energy losses due to the synchrotron (abbreviated as "SY" for 
subscripts appearing in the following equations) and IC are rather complicated, in addition 
to the energy dependent cross-section of the Klein-Nishina formula, 

/AE e \ I AE e \ I AE e 



where 



= e B (r), (12a) 



\ At 

1 lAE, 



w T El \ At /sy 

-^<^>«=£<£ww), < 12b > 

e i=0 

with w T = 1.018 x 1(T 16 cm 3 s _1 . See § 2.2 and Table 3 for [Tf, e^(r)] (i = 0-5), and A(E e , Tj) 
is given by equation (A6) in Appendix A, which comes from the Klein-Nishina cross-section. 
In Figure 4, we present the energy loss divided by E 2 e at the SS against E e separately for 
individual (virtual) photon fields as well as for superposed ones, — (AE e / At)f ic /E 2 , where 



we assume B± = 5/j,G, corresponding to e® = 0.93 eVcm 3 , for the magnetic field, and use 

e ph( r ©) P 

gas field. 



e ph( r o) presented in the second line from the bottom of Table 3 with r = r for the photon 
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Fig. 5.— Energy losses per unit time of CR electrons in ISRF and ISM at SS 
(r = r ) as a function of electron energy, shown separately for four processes, syn- 
chrotron +IC, bremsstrahlung, reacceleration, and ionization. We present three curves for 
synchrotron + IC, with the solid one from the Klein-Nishina cross-section, the broken one 
from the Thomson cross- sect ion, and the dotted one from the empirical one. 

For E e < 1 GeV, A(E e , 7$) 1, i.e., the Thomson cross-section is valid, so that equation 
(11) is separable in r and E e , leading to a simple expression, e(r)w T E^, with e(r) = e B (r) + 
Yli=o e ph( r )- 111 P rac tice, we find it is well reproduced by the following form over a wide 
energy range 

-(^) sic ^e(r)w T E 2 e - s ; 5 = 0.075, (13) 

while e(r) depends on E e very weakly. 

In Figure 5 we demonstrate the energy loss of individual processes separately, those 
due to "rad", "ion" and "sic" at the SS against the kinetic energy of the electron E e with 
e Q = 2 eVcm -3 and n = 1 H atoms cm -3 , where we plot the above empirical relationship (13) 
(dotted curve) and the energy gain due to the reacceleration (= "rea"; see next subsection) 
together. One finds that it reproduces satisfactorily the exact one (11) with equations (12a) 
and (12b). 

3.2. Energy gain due to the reacceleration 

In Paper II, we present the energy gain per unit time due to the reacceleration 



i(^L =n(r)wrca[Ee/GevrQ ' 



(14) 
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with 



Wrea = cCo; Co ~L^*' (15) 



where w rea , = 15.0 x 1CT 16 cm 3 s^ 1 in the case of, for instance, Co = 50 millibarn (mbarn), corre- 
sponding to the choice of a parameter set with v M = 20-30kms _1 (Alfven velocity), Uq = 0.06- 
0.14H atoms cm -3 , and Dq = 2x 10 28 cm 2 s _1 . The smallness of the gas density with <^ 1 
H atoms cm -3 indicates that the reacceleration process occurs even at some distance from 
the GP. 

The fluctuation in the energy gain due to the reacceleration is given (Gaisser 1990; 
Paper II) by 

1 iAEj\ 1 



Ef\ Ai A-,,, 2 



n(r)w rea [E e /GeV]- a . (16) 



3.3. Total energy loss and gain 

As discussed in the last two subsections, we have the total average energy-loss and the 
energy-gain per unit time 

"(^rL = < r ^ E ^ + c ( p W(30» ( 17 ) 

where 

W n (E e ) = w r&d (E e )E e + w ion (E e ) - w rea E l e ~ a , (18a) 

and 

W e (E e ) ~ w T E 2 ~ 5 - 5 = 0.075. (18b) 



One might note from equation (18a) that there exist two energies, E~ and E+, at which 
the first term proportional to n(r) in the right-hand of equation (17) becomes null. In Figure 
6, we demonstrate W n (E e ) against several choices of and find [E~, E+\ [0.1, 7] GeV in 
the case of Co = 50 mbarn. Namely, the synchrotron-IC is dominant for Eg^E^t, while the 
reacceleration is effective for E~ < E e < E^~, and the ionization for E e < E~. 

As discussed in § 2, the total number density of the ISM gas, n(r), and the total energy 
density of the ISRF, e(r), have complicated spatial distributions coming from local irregu- 
larities, which are not yet well established. On the other hand, in our previous papers, we 
have assumed a simple exponential- type form for n(r), smearing out the local irregularities, 



n(r) = n exp[-(r/r„ + \z\/z n )], 



(19) 
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Fig. 6. — Numerical values of W n (E e ) = w iad (E e )E e + w ioD (E e ) — w rca El a against E e for 
Co = 40-55 mbarn, where the screening effect for u> ra d is taken into account. 

greatly simplifying the complicated distributions given by equations (1), (2) with Tables 1, 
2, where no is the (interpolated) average gas density with approximately 1.5 H atoms cm -3 at 
the GC, and [r n , z n ] « [20, 0.2] kpc. 

In spite of such a simplification, we have found that our model reproduces remarkably 
well the experimental data on hadronic components. This tells us that charged CR com- 
ponents are well mixed during their propagation in the Galaxy over a residence time of 
approximately 10 7 yr, effectively smearing the local inhomogeneous structure of the ISM. In 
fact, it is well established that the anisotropy amplitude of CRs is of the level of at most 10~ 3 
at energies of 1-100 TeV (Sakakibara 1965; Nagashima et al. 1989; Cutler & Groom 1991). 
This is the reason why even the simplest leaky-box model and/or the simplified diffusion 
model such as, for instance, constant gas density and constant diffusion coefficient without 
spatial gradient, reproduces the CR hadronic components so well (Berezinskii et al. 1990). 

Now, corresponding to the simplification (19) for n(r), we assume the following simple 
exponential type form for e(r) as well 



where eo is the (interpolated) average energy density of the ISRF at the GC, and two param- 
eters, r e and z € , correspond to the scale heights for the spatial gradients, almost independent 
of the energy. Typically [e ; r e ,z e ] « [16eVcm~ 3 ; 4 kpc, 0.75 kpc] (Ishikawa 2010). 

However, while the simplifications given by equations (19) and (20) are applied for elec- 
trons (and hadrons), we stress here that those presented in §2 are actually used for D7's 
as discussed in §6, namely n h (r) with 11 h" = HI, HII, H 2 for n(r), and e^l(r) (2 = 0-5) for 




(20) 
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e(r) with weak energy dependences in as presented in Table 3. This is because D7's 
produced by CR hadrons and electrons are directly affected by the environment of ISM and 
ISRF around the birth site of the produced 7's. 



4. Diffusion equation for electron component 

4.1. Basic equation 

The transport equation for the electron density, N e (r; E e ,t), is given by (Berezinskii et 
al. 1990), 

^ - V • D(r; E e )V + A E j -N e (r; E e , t) = Q(r; E e , t), (21) 

with 

see equations (19) and (20) for the average energy-loss (-gain) in the all processes, with 
the replacement of [n(r),e(r)] in equation (17) by [n(r),e(r)], and equation (16) for the 
fluctuation of the energy gain in the reacceleration process respectively. For the diffusion 
coefficient and the source spectrum, we assume (note v ~ c, and R e ~ E e ) 

D(r;E e ) = E«D(r), Q(r;E e ,t) = E^Q(r;t), (23) 

with 

D(r) = D exp(r/r D + \z\/z D ), (24a) 

Q(r; t) = Q (t) exp[-(r/r Q + \z\/z Q )\. (246) 

In Table 4, we summarize parameters related to the scale heights, r D ,z D , . . . , which often 
appear in the present paper. 

Now, remembering W € (E e ) S> W n (E e ) in the high energy (HE) region, say, E e >E+ 
(^7GeV), and vice versa in the low energy (LE) region, Ee^E^t, the energy loss given by 
equation (17) is written as 



Age 



e{r)W e (E e ) + 0[n(r)W n {E e )]; E e > E+, (25a) 
n(r)W n (E e ) + 0[e(r)W e (E e )]; E e < E+, (25b) 



- 16 - 



Table 3. Summary of the numerical values of [Aq (/im), Tj(K), <7j], and those of 
e P *ho(eV cm ~ 3 ) for r > 3kpc in the second line from the bottom, and e^(eVcm -3 ) for 
r < 3kpc in the bottom line, where A 5 \ T 5 (K), and <x 2 have weak r-dependence as shown 
in remarks with r in kpc, while eg are independent of r except e^. The numerical value of 

(2) 

13.9 in e ph corresponds to the energy density of the stellar radiation for the population 
i = 2 at GC (see Figure 2). See also equation (7) for egg an d ^(r) ^ or r — 3kpc 



i = i = 1 i = 2 i = 3 i = 4 i = 5 remarks (r-dependence) 



Aq- 1 1.06(+3) 1.45(-1) 

Ti 2.73(+0) 9.92(+4) 

Gi 3.39(-l) 

s% 2.61(-1) 5.20(-l) 

eg 2.61(-1) 2.03(-l) 



1.05(+0) 7.50(+0) 2.00(+l) 

1.37(+4) 1.92(+3) 7.19(+2) 

6.10(-1) 2.20(-l) 2.19(-1) 

7.92(+0) 8.10(-1) 3.40(-l) 

1.39(+1) 3.17(-1) 1.33(-1) 



9.00(+l) 


Af(r) 


= 90e r / 37 - 5 


1.60(+2) 


r 5 (r) 


= 160e- r / 37 - 5 


4.72(-l) 


o- 2 (r) 


= 0.61e r / 87 - 


4.10(+0) 




_ (0 p -r/3.2 


1.61(+0) 


4 2 hW 


= 13.9e- r / 2 - 



Table 4. Summary of parameters often appearing in our propagation model, classifying 
them into two groups, one related to the gas density of ISM, n(r), and another to the 

energy density of ISRF, e(r). 



Parameters for ISM Typical values for ISM Parameters for ISRF Typical values for ISRF 



v = 



1 + z D /z n 



v = 0.04 - 0.06 



K = 



1 + z D /z e 



11/1 1 

2 \r D r n 



f n = [20 - 40]kpc 



-L = 1 ( — + — ) z n = [0.3 - 0.5]kpc 



2 y 2» 

1 1 



Zq 2z r 



z n w„ = 0.8-1.2 



1 



-(— - 

2 \r D + r £ 



11/1 1 

z e 2 \Zd z e 



k = 0.15 - 0.21 



U v = 2\fu + v 2 U u = 0.4 - 0.5 U K = 2\]k + k 2 U k = 0.8 - 1.2 



f e = [5 - 10]kpc 
z e = [1.0 - 1.4]kpc 
w K = 4.5 - 5.5 
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so that in the following discussion, we give first the solution of the diffusion equation (21) in 
the HE region, regarding nW n as a perturbative term, where we can neglect the fluctuation 
term due to the reacceleration. Next we give the solution in the LE region, regarding eW e 
as a perturbative term by contrast, which is completely the same as the former one after 
replacing [e, W e ] with [n, W„] (and vice versa), while we have to take the fluctuation term, 
(AE^) rca , into account in this case. 

Thus for the steady state (d/dt = 0), the solution of equation (21) in the HE region is 
devided into three 

N e , e ciN($ + N$+NM, (26a) 

where the first term is a principal one coming from eW e , the second term corresponds to the 
perturbative term from nW n , and the third term to the fluctuation due to the reacceleration 
given by the second term of the right-hand in equation (22), while it is negligible in practice, 

ivtV « 0. 

The solution in the LE region is similarly given by replacing the suffix e with n (and 
vice versa), but we can not neglect the fluctuation term Neln in contrast, 

N e>n c±NW+NM + N£l (26b) 
The first term in equation (26a) is written immediately as 

POO 

N^(r;E e ) = n^(r;y)f^(y;E e )dy, (27) 
Jo 



where ili - 1 and satisfy, 



e(r)c^-V-D(r)v]-nP{r-,y) = Q(r)5{y), (28a) 

cE: h ~ w^ {Ee) \ ' fi ° ){m Ee) = °' (28b) 



with/i 0) (0;£ e ) = E^- a . 



4.2. Solution in the steady state 

It is possible to solve exactly equation (28a) with use of the procedure presented in 
Paper I, after replacing n(r) by e(r), and we present here only the critical term related to 
(r, z; y), omitting constant terms such as Q and e (see Appendix B for the full form), 

IlW(r;y) oc exp[-s r y - \z\/z D ], (29) 
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+ k = t~ ( 3 °) 

e r cz^ \ k/ 1 + z D /z € 

with L> r = D(r, 0), e r = e(r, 0). As ili - 1 is of the form of e~ Sr ^, the Laplace transform of 
with respect to y, F$(E e ), is sufficient for our purpose to obtain the electron density, 

POO 

F^{Ee)= / e-^f^\y-E e )dy, 
Jo 

thus we have immediately from equation (28b) 

F ^ m = vue~7) J E dE ° E ^ e ~ ¥rAEe ' Eo) > ( 31 ) 

With f E £a 



Y„(E c ,E )=cs r J —dE. 



In the HE limit, E e 3> 1 GeV, using equation (18b), we find 

cE -{i+i-5) |- 1 + a _ 6 _ 2 



F^(E e )~ (1 _ £ a _ 6)w 



1 ~~ FT"/ + ••■ 



giving a spectral index with 7 + 1 — 5, where 5 (= 0.075) comes from the effect of the Klein- 
Nishina cross-section. Practically, however, it must be softer than the above index because 
of the exponential cutoff with e~ Ee ^ Ecut in the electron injection spectrum somewhere around 
20TeV (Reynolds & Keohane 1999; Hendrick & Reynolds 2001; Yamazaki et al. 2006). 

Now the principal term, NjP] , in equation (26a) for the electron density in the HE region, 
E e>E+, is given by 

4°)(r;E e )ocFW(E e )e-l^, (32) 

while the perturbative term, Ne°l, is obtained by the use of the iteration method as presented 
in Appendix Bl, giving A" e ( ,n /Ai? ~ 10% with the first iteration for E e > 1 GeV at SS as shown 
in Figure 24a. In practice, we perform only the first iteration, neglecting the second and 
higher iterations. Full form of A e e (r; E e ) is given by equation (B3). 

The numerical procedure in the LE region is similar to that in the HE region mentioned 
above by replacing the suffix "e" with "n" (and vice versa), while we have to take into 
account the third term in equation (26b), Aj,n, corresponding to the fluctuation. We find 
again that the perturbative term, N^e , is obtained by the use of the iteration method as 
presented in Appendix B2, giving AjJ/AS ~ 10% with the first iteration for E e < 10 GeV, 
as shown in Figure 24b. 
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Fig. 7. — Contributions of the fluctuation in the reacceleration with [a©, Co] = [180, 50] mbarn 
at SS for several sets of the free parameter 77* given by equation (B6). 



On the other hand, the numerical procedure in the fluctuation effect due to the reaccel- 
eration is a little bit cumbersome, which is presented in Appendix B3. We give an example of 
the ratio, Nell/[N^1 + Ne?e], at SS for the first iteration in Figure 7, where 77* is the effective 
ratio of the energy density to the gas density defined by equation (B6), approximately with 
2eV. One finds that it is significant around 0.3-1.5 GeV in the case of Co = 50 mbarn, boost- 
ing the solution without the fluctuation, Ne°n + Ne°) , by approximately 25%. So we perform 
only the first iteration also for Ne\l in the LE region, as the contribution coming from the 
second and higher iterations is at most of the magnitude of a few % or less (Ishikawa 2010). 
Full form of N e ^ n (r; E e ) is given by equation (B15). 

Finally, we give the electron density covering all energies so that it continues smoothly 
at the energy E c between the HE and LE regions at the SS (r = r ), with E c « in 
practice, but not always E c = E+, 



N e (r;E e ) 



N, 



e,0 



F r ,e( E e) -\z\/z D . for £ > E 

TP (TP \ ' e — " 

F r ,n( E e) -\z\/z D . f E < E 

TP (TP \ ' e — c ' 



(33a) 
(33b) 



with 



d F &te {E e 



dE e F Q , e (E c ) 



dE e F Q>n (E c ) 



see equations (B4) and (B13) for F re (E e ) and F r ^ n (E e ) respectively, and A^ e is determined 
by the normalization with the experimental data as discussed in § 5. 
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4.3. Numerical results 



In Figure 8 we show the numerical results of N e (r & ; E e )/N e $ in two cases, (a) [Co, 0©, 
s e }= [50, 180, 30 eV _1 ]mbarn with a = | (reacceleration with Kolmogorov-type spectrum 
in hydromagnetic turbulence), and (b) [0, 90, 15eV _1 ]mbarn with ct = \ (no reacceleration 
with Kraichnan-type spectrum) for e©/n = 2eV with /3 (= 7 + a) — 2.6, 2.7 and 2.8, see 
equation (30) with r = r Q for s , where we assume E cut = 20 TeV (Reynolds & Keohane 1999; 
Hendrick & Reynolds 2001; Yamazaki et al. 2006) in the electron injection spectrum with 
E^e~ Ee ^ Ecut , and the results show the use of both Klein-Nishina (solid curves) and Thomson 
(dotted curves) cross-sections. 

We find two critical points in Figure 8. First, those by the former cross-section give 
approximately 40-50% (20-30%) larger than those by the latter at 1 TeV (100 GeV), where 
the density is normalized at E e = 10 GeV, leading to significantly harder spectra than those 
with the Thomson cross-section, as expected. Similar results are also recently reported by 
Delahaye et al. (2010), while their main purpose is to study the nearby sources of electron 
and the positron excess problem as well, which are outside the range of the present paper. 

Second, the reacceleration effect is significant in the energy region less than 10 GeV as 
compared to the curves without the reacceleration process. Unfortunately, however, it is 
difficult to observe such a signal in the direct experimental data on the electron component 




Fig. 8. — Numerical results of the electron density, N e (r Q ; E e ), at the SS (r = r Q ) in the case 
of (a) reacceleration with [a; Co, Oq] = [|; 50mbarn, 180mbarn], and (b) no reacceleration 
with [|; 0, 90mbarn], for (3 = 2.6, 2.7, and 2.8, where the vertical axis is multiplied by E%, 
and normalized to E e = 10 GeV. We show results for two cross-sections in IC process from 
the Thomson (dotted curves) and Klein-Nishina (solid curves) formulae. 
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Galactocentric distance: r (kpc) Galactocentric distance: r (kpc) 



Fig. 9. — CR densities of (a) electrons, N e (r; E e ), and (b) protons, N p (r; E p ), as a function 
of the galactocentric distance r for several energies, where the vertical axis is normalized to 
the density at SS (r = r ). We present also those given by Hunter et al. (1997) (four kinds of 
square symbols), and by Strong et al. (1988) (thin filled histograms), where the CR densities 
are averaged azimuthally in each galactocentric quadrant by Hunter et al. (1997), assuming 
that they are coupled to the gas density of ISM independent of the energy, while azimuthal 
symmetric 7-ray emissivity is used by Strong et al. (1988). 

because of the modulation effect in the low energy region <j 5 GcV, which masks the electron 
flux boosted by the reacceleration, see Figure 13, even if it occurs actually. 

Next we examine the spatial dependence of the electron density. First it attenuates 
exponentially with the latitudinal distance z from the GP with the latitudinal scale height 
of the diffusion coefficient, ^ = 2-4 kpc, independent of the energy E e . This is the same 
result as in the case of the proton density, N p (r; E p ) (Paper V), namely the ratio of electron 
density to the proton density is independent of z. 

Contrary to the latitudinal behavior, the longitudinal behavior of the CR densities, 
N e (r; E e ) and N p (r;E p ), are somewhat complicated, both of which depend on the energy, 
and appear implicitly in the form of s r and a r (see Paper I for a r and its physical meaning). 
We present these in Figure 9 against the radial distance r for the (a) electron and (b) 
proton components, both normalized at the SS for four energies, 0.1, 1, 10, and 100 GeV, 
with (5 = 2.7, where the scale heights are set as [f n , f e ] = [30, 8] kpc and [z D ; z n , z e ] = [3; 0.2, 
0.75] kpc (see Table 4 for f n and f e ). We plot the results of Hunter et al. (1997; square symbols) 
and Strong et al. (1988; thin filled histogram) together, where the former are based on the 
assumption that the CR density is coupled to the density of ISM, and plotted separately for 
four galactocentric quadrants, I, II, III, and IV. We find that the radial dependence of the 
electron density, N e (r;E e ), is much stronger than that of the proton density, N p (r;E p ), in 
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the energy region of 1-100 GeV as expected, while the other two authors assume no spatial 
dependence in the energy spectrum, namely the shape of the energy spectrum at the SS is 
the same everywhere in the Galaxy. 



5. Electron-induced 7-ray spectrum 

For convenience in the following discussion, we summarize two cross-sections in Table 5, 
a EB (E e ,E y ) and a IC (E e , E y ; E ph ), each for the bremsstrahlung (abbreviated as "EB" for 
subscript attached here and in the following) and the IC processes respectively, where E p h is 
the energy of target photon before scattering. In these cross-sections, we take into account 
the screening effect for the bremsstrahlung (Koch & Motz 1959; Gould 1969), and the Klein- 
Nishina cross-section (Jones 1965, 1968; Blumenthal & Gould 1970) for IC. In the following 
discussion, we put E r (E e ) = F r>€ (E e ) for E e > E c , and E r (E e ) = F r ^ n (E e ) for E e < E c in 
equation (33) for simplicity. 

First we consider the emissivity of 7's from the bremsstrahlung at the position r, which 
is immediately written down as 



/•oo 

q EB (r;E 7 )= N e (r; E e )[n(r)ca EB (E e , E^)}dE e , (34) 
Je 1 

where the electron density, N e (r; E e ), is given by equation (33). For the numerical calculation 
of equation (34), we need the absolute electron density at r. To do so, we use the observational 
data on the electron intensity at the SS, dI®/dE e , which is related to the electron density 

by 



In practice, we normalize the electron density at E s = 10 GeV with use of the most 
recent data (see Fig. 14), where the solar modulation effect is negligible, 

cN® = cN e (r Q ;E s ) = 2.26m-V 1 GeV- 1 , 

corresponding to dIf/dE e = 0.180 m _2 sr _1 s _1 GeV -1 at E e = 10 GeV in Figure 14, while 
E s = 100 GeV (per nucleon) for the hadron-induced 7's (n° — > 27) with cN p (r Q ; E s ) = 
6.16m _2 s _1 GeV _1 (Paper V). One should keep in mind that the uncertainty in the normal- 
ization is of the magnitude as large as 10%. 



-23 - 



Table 5. Summary of the production cross-sections of 7-rays in the bremsstrahlung and 
the IC processes with x = E^/E e , where E e is the incident energy of electron, and E y is the 

energy of the produced 7's, and E ph the energy of the target photon before electron 
scattering. For the bremsstrahlung process, we present the cross-section in the case of only 
one-electron atoms (Z — 1), see Gould (1969) for two-electron atoms (Z — 2). 



bremsstrahlung (EB) 



inverse Compton (IC) 



2 2 \ 

a^=Aa f Z{Z + l)(-^) ; a f 



137 
2 



X\ dy 

y) y 



EB (x, E 7 ) = {l + (1 - x) 2 }Mx) - 3(1 - x)Mx) 
MX) = I + J\(y){l + 3^ (l + In ^) - 4 

My) = 1 - 



x 3 1 rfy 



y 3 J y 



[l + y 2 /(2a f Z) 



X = x(x,E y ) = 



1 1 

x m e c 
l-x 2E~, 



cr IC (E e , Ey, E ph )dE^ = 4>ic(x, q)^r 



(Tic 1 = 3<t t = 8ir 



2 \ 2 



b lc {x,q) = 2q\rvq+{l -q)[l + 2q + 



q = q(x,X) = 



1 x 2 



21 + x 



X = X(E e , E p h 



k = k(E ph ,T ph ) 



x 1 
1-xX 

k _ EphE e 



0| [m e c 2 /2] 2 



P h 



9 e = 6 e (E e ,T ph ) = 



kvTph 
m e c 2 /2 



y/kBTphEe 



Thus taking care of the terms related to r, we have 



q w (r; E 7 ) 
n(r)wgiN° 



1 



F r (E x ) dx 
F Q (E s )x^ 



EB (a;,-E 7 ) 



(35) 



for E 7 > E c with ^ = i^/a, and w^l — ca^l = 1.39xl0 _16 cm 3 s _1 for the hydrogen gas 
(Z = 1), where one should take care of the energy range E e < E c in the case of E 1 < E c . 

Next we consider the emissivity of 7's coming from the IC process, which is somewhat 
complicated, as there are several kinds of target photons with different energy density as 
well as with different scale heights in the spatial gradient. Here we present a result only, 
taking into account the six wavelength bands in eg(r) (i = 0-5) (see eq. [7] and Table 3), 



where $ic (x,E 7 ) is given by equation (C2), see Appendix C for the details. 

Let us show the numerical results for two cases of emissivity in Figure 10, (a) r/r = 
0.5,1,2 with z = in the GP, and (b) z = 0.2, 0.4, 0.6kpc with r = r normal to the 
GP at SS, assuming (3 = 7 + a =2.7, where we present separately those coming from tt° 
(solid curves), EB (broken curves), and IC (dotted curves). One finds that EB-7's and 7r°-7's 
are comparable around 50 MeV, and IC-7's and 7r°-7's around two energies, ~ 20 MeV and 



See Paper V for the emissivity originating in tt°, q w o(r; E 7 ), while we use more realistic 
gas density, n(r), in the present paper. Note also in q n o that the semi-empirical production 
cross-section of 7's, cr pp ^ 1 (E p , E^), in proton-proton collision we use is valid over very wide 
energy ranges, 1 GeV-1 PeV, reproducing nicely various kinds of physical quantities such as 
psuedo-rapidity, energy spectrum, multiplicity, etc, obtained by both the accelerator and CR 
experiments with local target layer (Suzuki, Watanabe & Shibata 2005). 

Once we have the emissivity of 7's induced by the interaction between the electrons and 
the media of ISM and ISRF, we can obtain immediately the intensity of 7's observed at the 
SS (r = r ), coming from the direction 6(1, b) 




(36) 



ITeV. 



d^l^O-E,) _ 1 r 

dEjdld(sin b) An J 




with 



5 



g 7 (r; f? 7 ) = q EB (r; £ 7 ) + ^ q$(r; E 7 ), 
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10 1 10° 10 1 

y-ray energy; £ r (GeV) 



10 1 10° 10 1 10 2 

y-ray energy; E r (GeY) 



Fig. 10. — Emissivity of 7's (a) at three radial distances, r/r Q = 0.5, 1.0, 1.5 in the galactic 
plane (z — 0), and (b) at three latitudinal distances, z = 0.2, 0.4, 0.6 kpc with r = r , where 
the three components for 7's emission, 7r° (solid curves), EB (broken curves), and IC (dotted 
curves) are shown separately. 



where the integration with respect to s is performed along the arrival direction of 7's, 6(1, b), 
at the SS, and r(r, z) is bound to (s; I, b) as follows, 

r(s; l,b) = \Jr£ + s 2 cos 2 b — 2r Q s cos b cos I, 
z(s; b) = ssinfr. 



6. Comparison with the observational data 

6.1. Critical parameters 

We assume that the source distribution of electron component, Q(r;E e ), is the same 
as that of the hadronic component except for the cutoff electron energy, for instance E cnt 
20 TeV, with the supernova remnants as the main energy supply, while the pulsars and pulsar 
wind nebulae might contribute to them as well, particularly to positrons and electrons (for 
instance, Delahaye et al. 2010). So the galactic parameters used in the present work are 
essentially the same as those appearing in Papers I-V, and we summarize them briefly in 
the following. 

The recent observational data on the energy spactra of CR hadronic components give 
indices with 2.74 ± 0.08 for proton (Derbina et al. 2005), and with a common value of ~2.7 
for nuclei between the oxygen and iron (Miiller 2009), whereas there still remains uncertainty 
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for helium, for instance with 2.68 ± 0.05 by JACEE (Asakimori et al. 1998) in contrast to 
2.78 ± 0.20 by RUNJOB (Derbina et al. 2005). Note that PAMELA (Picozza et al. 2007) 
reports recently a common index of 2.73 in both the proton and helium spectra, albeit the 
energy region is limited below 500 GeV. Any way, the spectrum index (3 of proton, must lie 
well within 2.7-2.8 in the high energy region at the SS, which is the most effective element 
for the hadron-induced D7's. See Paper V for the contribution of helium and nuclei to 
D7's, which is taken into account by introducing the enhancement factor with 1.53. So in 
the present paper we use the critical parameter (3 in place of 7 (source index of the energy 
spectrum) with f3 = 7 + a, and consider three values of (3; 2.6, 2.7, and 2.8, each for a — | 
(Kolmogoroph-type spectrum) and \ (Kraichnan-type spectrum). 

There are three galactic parameters, [D(r), n(r), Q(r)], in our approach to the CR 
propagation, and six scale heights for longitudinal and latitudibal directions correspoding to 
each one, [r D ,r n ,r Q ] and [z D ,z n ,z Q ], respectively. In practice, however, explicit parameters 
needed to compare with the experimental data appear in two critical ones alone, a r and Co, 
besides [a, f$\ mentioned above, while the parameter, \x r = 2z D / y/f D r , is also important for 
the study of the CR isotopes (f : normalized life time of an isotope with 10 6 yr). 

For electron components, the additional parameter newly appears, s r , given by equation 
(30), physical meaning of which is essentially the same as a r ; i.e., while the inverse of a r 
gives the average path length, x r , in units of cm -2 in ISM as discussed in Paper I, that of 
s r corresponds to the average path length, y r , in units of eVcm~ 2 in ISRF, namely the total 
amount of photon- gas energy that CR has passed through the ISRF. 

Now from equation (33), one should remark that there appear only three critical pa- 
rameters, [(0, <Tq, Sq], in F r>€ (E e ) and F r ^ n (E e ) needed to compare with the observational 
data, aside from two critical indices, [a, 0\, note that various galactic parameters such as 
the diffusion constant, gas density, energy density, their scale heights, etc are all involved 
implicitly in these three ones. 



6.2. Charged components 

6.2.1. Hadron components 

As we have presented the experimental results on CR hadron components in the past 
papers (Papers I-IV), we give here only three kinds of secondary-to-primary ratio with new 
data, B/C, sub-Fe/Fe, and p/p, that have since become available. See Paper III for the 
secondary unstable nuclei, while new data are still not available. 
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In Figure 11, we present B/C and sub-Fe/Fe, plotted together with new ones from 
CREAM (Ahn et al. 2008) and TRACER (Miiller 2009), where we plot also RUN JOB 
(Derbina et al. 2005) data for reference, while the data quality is rather poor with large 




kinetic energy : £ N [GeV/nucleon] kinetic energy : E N [GeV/nucleon] 

Fig. 11.— Energy dependence of the secondary-to-primary ratio for boron/carbon and 
sub-iron/iron. See Paper II and references therein for the experimental data, while CREAM 
(Ahn et al. 2008) and TRACER (Miiller 2009) data are newly plotted. Numerical curves 
are demonstrated for two cases; (a) reacceleration with (a, £q) = (|, 50mbarn) and (b) no 
reacceleration with (~, 0). 




kinetic energy : E F (GeV) 



Fig. 12.— Energy dependence of the antiproton-to-proton ratio, where we assume the 
reacceleration model with a = | for six modulation parameters, 0.2, 0.3, 0.5, 0.7, 1.0, and 
1.5 GV, and A < (A > 0) corresponds to positive (negative) polarity state in heliospheric 
magnetic field, although the present calculations do not take the effect into account. See 
Paper IV and the references therein for the data other than PAMELA (Adriani et al. 2010). 
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Fig. 13.— Electron energy spectra in two models, (a) with reacceleration and (b) no 
reacceleration, compared with the measurements, where the vertical axis is multiplied by 
E^. All numerical values are normalized to _E e = 10GeV, and several sets of (/3, s r ) are 
assumed. See text for references for individual experimental data. 

atmospheric correction. We compare our numerical results with the data for two models, 
(a) reacceleration with the set of [Co, 0"©] = [50, 150-300] mbarn, for a = |, and (b) no reac- 
celeration with [0, 75-150] mbarn, for ci = \- It is still not clear which model reproduces 
the experimental data more satisfactorily. As is well known, the advantage of the former 
explains naturally the drop of the ratio in the lower energy region around ACE/CRIS (Davis 
et al. 2000) without assuming an ad hoc drop in the path length distribution. 

Next we present p/p in Figure 12, plotted together with new data from PAMELA 
(Adriani et al. 2010), where we present numerical curves with several sets of modulation 
parameters, 0.2-1.5 GV, for the reacceleration model shown in Figure 11a. One finds that 
our result is in good agreement with the PAMELA in the high energy region around 100 GeV, 
where the modulation effect is absolutely negligible. 



6.2.2. Electron component 

Let us present the electron data separately before and after FERMI, where "electron" 
denotes both electron and positron. First in Figure 13 we present the electron energy spec- 
trum before FERMI, where the experimental data are presented for those reported in the 
period from 1994 to 2008 alone (Golden et al. 1994; Kobayashi et al. 1999; DuVernois et 
al. 2001; Torii et al. 2001, 2006; Aguilar et al. 2002; Chang et al. 2008), and also plotted 
are the data {filled purple squares) for reference after applying a demodulated correction to 
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HEAT data, HEAT-LIS, (DuVernois et al. 2001) using the force-field approximation with 
the modulation parameter of 755 MV (670 MV) for the 1994 (1995) data. 

The numerical curves are normalized at 10 GeV with two indices, = 2.7, 2.8, assuming 
two models, (a) reacceleration and (b) no reacceleration each with the same parameter sets 
as those used in Figure 11, while we assume additionally two cases of s Q , [20, 30] mbarn for 
the reacceleration (a), and [10, 15] mbarn for no reacceleration (b). Aside from the prominent 
spectral features around 500 GeV appearing in ATIC (Chang et al. 2008) and PPB-BETS 
(Torii et al. 2006) data, our model with the reacceleration reproduces the data well in the 
higher energy region, >10 GeV, in Figure 13a, where the solar modulation effect is small. 
On the other hand, the model without reacceleration in Figure 13b is somewhat difficult to 
fit to the demodulated HEAT-LIS data. 

Now, in Figure 14 we present the most recent data obtained by FERMI (Abdo et al. 
2009, 2010b) and H.E.S.S. (Aharonian et al. 2009) together with those presented in Figure 13, 
where numerical curves are the same as shown in Figure 13a. We find that both FERMI and 
H.E.S.S. data do not exhibit the prominent bump around 500 GeV reported by ATIC and 
PPB-BETS, with both giving a spectrum falling with energy as E~ 3 up to ITeV, which is 
not inconsistent with emulsion chamber data (Kobayashi et al. 1999) within the statistical 
errors. Looking Figure 14, however, we find that FERMI and H.E.S.S. data seem to deviate 
systematically from numerical curves with an enhancement by 20-30% around 500 GeV, 
indicating still some additional local sources of high energy CR electrons, which will be 
discussed again in § 7. 




Fig. 14.— Same as Figure 13a, but with FERMI (Abdo et al. 2009, 2010b) and H.E.S.S. 
(Aharonian et al. 2009), where drawn are numerical curves with the reacceleration shown in 
Figure 13a. 
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Fig. 15. — An example of the estimation of the BG [black horizontal line) using the latitu- 
dinal D7's data from EGRET (Hunter et al. 1997) with the energy interval of 300-500 MeV 
averaged over the whole radial direction, Z = 0°-360°. 

6.3. Diffuse 7-ray component 

6.3.1. Isotropic background 7 -rays 

D7's near the GP are mainly hadron-induced 27) and electron- induced (EB + IC). 

In addition to these two components, we have isotropic background 7's (BGs) with various 
origins such as extragalactic sources (EGs), unidentified sources, instrumental sources, dark 
matter (DM), etc, so that the BGs depend on individual detectors with different sensitivity 
in energy and the angular resolution, while depending on the propagation model as well. 
Therefore it is not easy task to estimate the extragalactic D7, while its origin is one of the 
fundamental problems in astrophysics, studied in so many papers with various candidates; 
unresolved blazers (e.g. Stecker & Salamon 1996; Chiang & Mukherjee 1998; Miicke & Pohl 
2000), intergalactic shocks produced by the assembly of large-scale structures (e.g. Loeb & 
Waxman 2000; Totani & Kitayama 2000; Miniati et al. 2000; Gabici & Blasi 2003), dark 
matter annihilation (e.g. Bergstrom 2000; Ullio et al. 2002; Ahn et al. 2007), etc. In the 
present paper, however, we use the acronym "BGs" all together for D7's other than those 
induced by ir°, EB (bremsstrahlung) and IC, while acknowledging EGRET and FERMI 
teams have estimated very carefully the EG-7 intensity. 

In Figure 15, we present an example of EGRET data (histogram; source subtracted) 
(Hunter et al. 1997) together with numerical curves on the latitudinal distribution aver- 
aged over full longitude ranges, 0°-360° with the energy interval 300-500 MeV, where we 
give the contributions of D7's separately from 7r° (solid red), EB (dotted red), IC (bro- 
ken red), BG (solid black), and total flux, 7r°+EB+IC+BG (heavy solid red), assuming 



-31 - 




Fig. 16.— The BG spectrum obtained by COMPTEL (Kappadath et al. 1996), EGRET 
(Sreekumar et al. 1998), EGRET (revised; Strong et al. 2004), and FERMI (Abdo et al. 
2010a), where also plotted are those estimated by the present work using EGRET and 
FERMI data (see Figs. 15 and 17). Dotted line is given by Abdo et al., and the solid curve 
is used in the present work (see text), modifying it slightly in the low energy region. 



[Co? s Q ] = [50, 180, 30 eV _1 ]mbarn with [a,/3] = [|, 2.7]. Here we draw a horizontal line 
for BG by the use of the least square method so that the histogram is well reproduced, 
where the fitting is applied for |6| < 60° as there remain considerable uncertainties in the 
latitudinal distribution for both the ISM and ISRF far distant from the GP, see E h (r,z) 
("h" = HI, H 2 ) in equation (2) and the scale height z p h in equation (7). It is remarkable 
that the numerical curve is in good agreement with the data not only in shape, but also in 
absolute value, except the high latitude around the galactic pole. 

In Figure 16, we summarize the intensity of BGs obtained by past works, Kappadath 
et al. (1996) for COMPTEL, Sreekumar et al. (1998) for EGRET, Strong et al. (2004) 
for EGRET (revised), and Abdo et al. (2010a) for FERMI, where also plotted are those 
estimated in this work (see Figs. 15 and 17) for the reference, six points (open circles) for 
EGRET and one point (filled circle) for FERMI. We draw a broken line given by Abdo et al. 
with dJ BG /d£7 7 = 9.6 • 10~ 3 x E~ 2Al in units of [cm-V^^MeV -1 ] with £ 7 in MeV, and 
a solid curve with dI* G /dE 1 = dI BG /dE 1 x [1 — exp (— 0.33-E^' 4 )], used in the present work, 
slightly modifying the FERMI result in the low energy region, while the modification does 
not affect any change for the results. Significant difference between EGRET and FERMI, 
with the former giving much harder spectrum than the latter, might be due to a diffrent 
model in CR propagation as well as those in ISM and ISRF. 



-32 - 




Fig. 17. — Latitudinal distributions of D7's obtained by a) EGRET (Hunter et al. 1997) 
with the energy range 300-500 MeV, and b) FERMI (Porter 2009) with the energy range 
350-480 MeV, both averaged over the the whole radial direction, / = 0°-360°, where drawn 
are three curves with /3 = 2.6 (green), 2.7 (red), and 2.8 (blue), taking the BG contributions 
into account. 



6.3.2. Spatial distribution 

We present two examples of the latitudinal distributions for EGRET and FERMI (Porter 
2009) with the energy interval around [300-500] MeV in Figures 17a and 17b respectively, 
together with our numerical results taking the BG contribution (broken- dotted lines) into 
account mentioned above, dI* G /dE y , where plotted are three curves for each figure with 
(3 = 2.6 (green), 2.7 (red), and 2.8 (blue). One finds the agreement between the data and 
the curves is excellent except the high latitude |6| >,60°. In these calculations, we take the 
angular resolution (PSF) effect with the energy dependence into account, for instance, with 
7° (HWHM) at 30-50 MeV (Hunter et al. 1997). 

Corresponding to the latitudinal distributions as shown in Figures 17a and 17b, we 
demonstrate the longitudinal distributions near the GP in Figures 18a and 18b, where nu- 
merical curves are shifted by Al = +10° in both EGRET and FERMI so that experimental 
data are reproduced more satisfactorily. Again we find the numerical results are in nice 
agreement with the data in both shape and absolute value, and consistent with /3 ~ 2.7. 
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Fig. 18. — Same as Figure 17, but for longitudinal distribution with the same condition, 
where BG contributions are not presented. 




j'-ray energy; E y (GeV) y-ray energy; E y (GeV) 

Fig. 19. — Differential energy spectra of D7's averaged over the whole radial direction with 
\b\ < 10° obtained by COMPTEL (Kappadath et al. 1996), EGRET (Hunter et al. 1997) and 
FERMI (Abdo et al. 2010b). Numerical curves are demonstrated for two ) with the 

reacceleration and b) no reacceleration, each presented separately for individual components. 



6.3.3. Energy spectrum 

First, in Figure 19a we present the energy spectrum of D7's averaged over the field of 
view with —60° < I < 60° and < 10°, where numerical curves with the reacceleration are 
also presented separately for those coming from it , EB, and IC (all with colored thin solid 
curves), BG (heavy black solid curve), and total (colored heavy solid curves). Here and in 
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Fig. 20. — Same as Figure 19a, but those averaged over three different ranges in the galactic 
latitude, a) low latitude with 10° < \b\ < 20°, b) mid latitude with 20° < \b\ < 60°, and c) 
high latitude with 60° < |6| < 90°. Parameter sets for the numerical calculations are the 
same as those used in Figure 19a with the reacceleration. 



the following we omit EGRET data in GeV region, because of the instrumental problem in 
detection of 7's (Stecker et al. 2008). One finds that the curve with (3 = 2.7 (red) is in good 
agreement with the data in the energy region below 1 GeV, while they deviate slightly from 
the curve above 1 GeV, with approximately 20% enhancement. 

Figure 19b reproduces Figure 19a, but for curves without reacceleration, corresponding 
to Figures lib and 13b (see also Fig. 8b). The fit is not as good as for the reacceleration 
model, particularly in the low energy region, E 1 ^ 200 MeV, with ~ 40% enhancement, while 
with ~20% in the high energy region, E 1 ^> lGeV, giving nearly the same enhancement as 
in the case of (a) with the reacceleration. 

Second, in Figures 20 and 21, we present the energy spectra of D7's for different sky 



views (Abdo et al. 2010a; see also supplementary material at http:/ /link.aps.org/supplement-al/10.1103/Ph 
Figure 20 shows those averaged over independent galactic latitude ranges covering low, mid 
and high galactic latitudes, a) 10° < \b\ < 20°, b) 20° < \b\ < 60° and c) 60° < \b\ < 90° respec- 
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Fig. 21. — Same as Figure 20, but those averaged over different hemispheres on the sky for 
the galactic latitudes with \b\ > 10°, which are centered at a) the north galactic pole {open 
square) and the south galactic pole {filled square), b) the galactic center with 270° < I < 90° 
(inner Galaxy), and c) the anticenter 90° < / < 270° (outer Galaxy). 

tively. Figure 21 shows those averaged over different hemishperes, which are, a) centered 
at the north {b > 0°; open squares) and south {b < 0°; filled squares) galactic poles, b) the 
galactic center (270° < / < 90°), and c) anticenter (90° < I < 270°), all with the galactic 
latitudes excluding |6| < 10°. In these figures, we subtract 7's coming from point sources 
based on the FERMI catalog. 

It is remarkable in Figure 20 that EGRET and FERMI data agree pretty well with 
each other, overlapping nicely around 0.2-lGeV, in all latitude ranges. One finds that the 
numerical curves with /3 = 2.7 reproduce generally well both the EGRET and FERMI data in 
Figures 20 and 21 but 20c, taking account of the uncertainties in various galactic parameters, 
particularly in those related to the ISM and ISRF. 

On the other hand, in Figure 20c for the high latitude, |6| > 60°, we have a noticeable 
enhancement in FERMI with approximately 70% as compared to the numeical curves. To 
see the deviation more clearly, we present them all together in Figure 22, where we show ad- 
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Fig. 22. — Same as Figure 20, but two cases of numerical curves, a) with the BG-7 given by 
FERMI (Abdo et al. 2010a) (solid curves), and b) with the BG-7 by EGRET (Sreekumar 
et al. 1998) (dotted curves). 




y-ray energy; £ r (GeV) 



Fig. 23. — Same as Figure 21, but two cases of numerical curves, a) with the BG-7 given by 
FERMI (Abdo et al. 2010a) (solid curves), and b) with the BG-7 by EGRET (Sreekumar 
et al. 1998) (dotted curves). 

ditionally numerical curves (dotted colors) using EGRET-BG obtained by Sreekumar (1998) 
for reference (see Fig. 16). Figure 23 reproduces Figure 21 with numerical curves using 
EGRET-BG (dotted colors) in addition to those using FERMI-BG (solid colors). 

One finds the spectrum shapes with EGRET-BG are quite different from the data in 
the high energy region, although the enhancement is rather improved in the energy region 
< 1 GeV, which is discussed again in the next section. 
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7. Discussion and summary 

We have studied the diffusion-halo model with stochastic reacceleration, comparing 
it with the most recent data on hadronic, electronic and D7 components. We have two 
particular interests: to find an unified model for the CR acceleration and propagation from 
the viewpoint of astrophysics, and to search for a signal of novel sources such as PBH 
and/or DM from the viewpoints of particle physics and cosmology. Both are of course 
closely connected with each other in the sense that the knowledge of the former is decisive 
in confirming the latter. While several groups (Torii et al. 2006; Chang et al. 2008) have 
reported the possibility of annihilation and/or decay of DM particles, giving a significant 
bump in electron flux around 500 GeV, FERMI (Abdo et al. 2009) and H.E.S.S. (Aharonian 
et al. 2009) give a rather flat spectrum up to 1 TeV without the prominent excess. In the 
present paper, however, we have focussed our interest rather conservatively on the internal 
consistency among various CR components from the view point of astrophysics, leaving the 
puzzle of the possible electron/positron-excess to further observations and mutual cross- 
checks in data analysis among individual groups. 

In our past works on the hadronic component, we concluded that the diffusion- halo 
model with the reacceleration with the parameter set, [Co, <7©] = [50, 180] mbarn with [a, f3] = 
[|, 2.7-2.8], is in harmony with the CR hadron data presently available. The most recent 
data on the B/C ratio by CREAM and TRACER (Fig. 11) as well as on the p/p ratio by 
PAMELA (Fig. 12) also support the present model. However, it is worth mentioning here 
that our interpretation for the energy dependence of the B/C ratio is somewhat different 
from that by CREAM (Ahn et al. 2008) and TRACER (Muller 2009). 

They claim that the index a favors 0.5-0.6 instead of |, resulting in a rapid decrease with 
energy for the interstellar propagation path length. In contrast to their interpretation, we 
would like to point out that the value of 0.5-0.6 is not fundamental, but is rather accidental 
due to the reacceleration effect, namely it is boosted upward around the GeV region by 
the energy gain, resulting coincidentally in the soft slope with 0.5-0.6 in the energy region 
1-100 GeV. The intrinsic one must be | (Kolmogorov-type for wave number spectrum in 
hydromagnetic turbulence), leading to 1) a natural drop in path length distribution in the low 
energy region < 1 GeV without introducing an artificial break there, as originally proposed 
by Simon et al. (1986), and 2) a reasonable amplitude in the anisotropy of CR's with the 
level of 10 -3 in TeV region nowadays established experimentally. 

We apply the diffusion-halo model with and without the stochastic reacceleration for the 
electron and D7 components. Apart from the electron-anomaly around 500 GeV, we find that 
the parameter set with the reacceleration, [a, (3;(o, <7 Q ], expected from hadron component 
reproduces rather well both the spectrum shape and the absolute value in both the electron 
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(Fig. 13a) and D7 (Figs. 17-21) components, assuming the additional parameter s with 
20-30 eV _1 mbarn. Physical meanings of the numerical set with [£ , <r ] = [50, 180] mbarn 
are discussed in Papers I III in connection with the diffusion constant D Q , gas density n Q , 
their scale heights, z D ,z n , etc, giving reasonable values matched with the observtional data. 

Let us consider the physical meaning of 20-30 eV _1 mbarn in s . The relation between 
<j and s Q is given by 

s _ n Q 1 + 1/k _ n Q 2 + z D /z e 
l + l/i/ e Q 2 + z D /z n ' 

see §4.2 for s and <r with r = r , and Table 4 for v and /c. Namely, it is closely related to 
the ratio of the energy density e to the gas density n at the SS, for the smeared energy 
density, smeared gas density respectively, and three latitudinal scale heights, z e , z n and 
z D . As discussed in §§2 and 3.1, we have n = n® + + n n = 1.14H atoms cm -3 , and 
e © = 6 b + e ph ~ 2.3eVcm -3 , leading to e /n m 2eV. Remembering that the scale heights 
used in the present paper are [z n , z e ; z D ] = [0.2, 0.75; 3.0] kpc, we find s = 32eV _1 mbarn for 
<7 = 180 mbarn, giving a consistent result, while the latter with 180 mbarn is expected from 
the relation, a e ~ D Q /[n e cz D z n ] with a reasonable set [D Q , n ] = [3x 10 28 cm 2 s _1 , 1cm -3 ] 
and [z D) z n ] = [3, 0.2] kpc as discussed in Papers I, II. 

As mentioned above, the electron spectra currently available are generally in agreement 
with those expected from the hadron spectra, considering the uncertainties inherent in both 
the experimental data and the numerical parameters, but not quite satisfactory, with the 
FERMI data giving the excess by 20-30% around several hundred GeV compared to the 
numerical results as seen in Figure 14. It might be related to the positron excess around 10- 
100 GeV observed by PAMELA (Adriani et al. 2009), indicating some nearby sources and/or 
exotic ones from DM annihilation or decay, while beyond the subject of the present paper. In 
fact most recently Delahaye et al. (2010) show that the electron spectra with FERMI, HESS 
and PAMELA are reproduced rather well by the standard astrophysical processes, assuming 
two sources separately, the distant and local nearby ones, whereas they stress that there 
remain too large theoretical uncertainties to build a standard model for CR electrons. So it 
is critical to study the D7's and diffuse radio emissions simultaneously in order to reduce the 
uncertainties inherent in the galactic parameters assumed for the numerical calculations. 

We compared our numerical results on the energy spectrum of D7's with EGRET and 
FERMI data for several sets of the field of view (Figs. 17-21), and found that overall, 
the CR data, hadron (Fig. 6 in Paper V) and electron (Fig. 14) components, reproduce 
rather satisfactorily D7's for both EGRET and FERMI, considering the fact that we have 
uncertainties with at least 10-20% in the galactic parameters assumed here as well as in the 
flux normalization of the hadron and electron components. Small enhancements of D7's in 
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GeV region (Figs. 19, 20), albeit they are still within the uncertainties, may indicate those 
from nearby sources such as the supernova remnants, pulsars, and pusar wind nebulae. 

We found, however, that FERMI data give the significant excess with approximately 
70% or more in the high latitude (Fig. 22), well beyond the uncertaities, against the numerical 
results in GeV region. This result may indicate a signature of very large electron-halo far 
distant from the GP, with, for instance, as large as 25kpc (Keshet et al. 2004), and/or 
something else coming from the cosmological origin. We are also concerned if the exess here 
discussed relates to those appearing in the electron spectrum between 100 and lOOOGeV 
observed by FERMI and HESS (Fig. 14) and in the positron spectrum around several tens 
GeV by PAMELA. To make clear the correlation between these excesses, D7's in the high 
latitude and the electrons/positrons around several tens to hundred GeV, crucially important 
is the anisotropy study for the high energy electron, which will be discussed elsewhere in the 
near future. 

Finally we briefly argue the electron spectrum obtained by FERMI from the observa- 
tional point of view, aside from the prominent bumps indicated by ATIC and PPB-BETS. 
FERMI is indeed excellent in the observation for 7-rays, we have some concerns about the 
separation of electrons from hadrons as well as their energy determination in the high energy 
region, while acknowledging the team have studied very carefully the reliability from various 
kinds of checks, with both beam tests and the simulational analyses. 

Nevertheless, one should keep in mind that FERMI is not purely- direct observations for 
electrons, but quasi- direct ones in the sense that electron events are selected by statistical 
analysis based on simulations for the spread of electron showers, where a small number of 
electrons are statistically selected from a large proton background. In contrast to these quasi- 
direct experiments, the PAMELA apparatus consists of a permanent magnetic spectrometer 
with a silicon tracking system, providing good identification between electrons and positrons, 
though limited to a maximum detectable rigidity (MDR) of 100 GV. 

Anyway, we await further studies and mutual cross-checks among the groups from var- 
ious points of view to get a firm conclusion for the electron-excess around 500GeV, while 
not so prominent as given by ATIC and PPB-BETS. So results from the AMS program 
(Bindi 2009), AMS-02, will be of particular interest. This program aims at high precision 
measurements of CR (both electron and hadron) and 7-ray fluxes from a few hundred MeV 
to a few TeV using a super- conducting magnet^, with the space shuttle launch scheduled for 
September 2010. We also look forward to D7's data from ground-based telescopes currently 



1 After submitting the present paper, we find that they decided to use the permanent magnet in place of 
the super-conducting magnet (Kounine 2010). 
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operating such as H.E.S.S., MAGIC, and VERITAS, as well as the CTA-program now under 
consideration (Caballero et al. 2008), the threshold energies of which are now overlapping 
with the FERMI satellite data. 

We are very grateful to P.G. Edwards (CSIRO Astronomy and Space Science) for his 
careful reading of the manuscript and valuable comments. 

APPENDIX A 
ENERGY LOSS OF ELECTRONS IN ISM and ISRF 

The energy-loss rate due to the bremsstrahlung in the gas density n(r) is given by 
/ AE \ f Ee 

\^T/rad = J E 1 [n(r)ca- E n{E e ,E 1 j\dE 1 = n(r)w TSL(i (E e )E e , (Al) 

w rad (E e ) = u4b J <Peb(x, E e x)dx, 

with ° 

W W = ca { ^ = 4ca f Z(Z + 1) = 1.39 x lO^cmV 1 , 

\m e c 2 / 

for hydrogen atoms, and see the left-hand side of Table 5 for (p EB (x, E y ). 

For E e 3> m e c 2 , we can use the complete screening cross-section, leading to the well- 
known result 

"(lI'L = n ( r ) w ™ E <> w ^ = 7 - 30 x io-^s- 1 . 

On the other hand the energy-loss rate due to the IC is given, taking into account the 
energy spectrum of the target photon at r, n p h(r; E p t), by 



where 



\ At he 

with 



/ AE \ f°° f Eyi 

1 -) = / dE ph / E^[n^{r\E^)ca lo {E^E^E^)\dE^ (A2) 



X 4E ph E e 

E M = E M (E e: E ph ) = E e — — -; X = X(E e , E ph ) = — , (A3) 

1 + X ' [m e c z Y 

see the right-hand side of Table 5 for a IG (E e , E y ; E ph ). Here we omit the suffix % introduced 
in §2.2 for simplicity. For E e ^> m e c 2 , equivalently X >> 1, one finds a reasonable result, 
E M m E e , leading to E ph < E y < E e . 
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From equation (4) in the text 

"eph(r) 



nph(»*;Sph)d-Eph 



"^B^ph 



where W ph (k) is the Planck function for the 2.7 K CMB, and the gaussian function given by 
equation (5) for the stellar radiation and the re-emission from the dust grains. 

The integration with respect to E 1 is given (Jones 1965, 1968) by, (see eq. [A3] for X) 

(0) 2 

E y a lc (E e , E y ; E ph )dE 7 ps a ic E e S IC (X)/X , 



L 



with 

a$ = 3cr T = 2.00 x 10~ 24 cm 2 , (cr T : Thomson cross-section), 

and 

'X 6\ W1 llX 3 /12 + 6X 2 + 9X + 4 /- x ln(l + t) 



o / v n / x „ 6\, >. 11X712 + 6X 2 + 9X + 4 n n r 
S IC (X)=(- +6+ -) H l + X) 2-2 1 



t 



dt. 



One should note that the approximation used above is only E e ^> m e c 2 , readily satisfying 
the condition in the energy region of interest, E e :> 10 MeV. 

Now, we have the energy-loss rate due to the IC scattering in a compact form after 
integrating over the energy E ph of the target photon in equation (A2), 

-(^f) IC = ^(r)w T A(E e ,T ph )El (A4) 
where E e is in units of GeV and e p h in eVcm -3 , and 

w T = ^Zrw = 1-018 x lO-^cmV 1 , (A5) 
3 [m e c 2 / GeV\ 2 

POD 

A(E e ,T ph ) = A(6 e ) = 9 / S lc (X)W ph (9 2 e X)dX/X A , (A6) 

Jo 

with 

6 e =0 e (E e ,T ph )= (AT) 

The above discussions are applicable also for the synchrotron radiation, since it is caused 
by the collision between an electron and the virtual photon induced by the magnetic field. 
Practically, however, we have the condition hw c r e <C m e c 2 (r e : Lorentz factor of elec- 
tron) with u c = eHj_/m e c 2 , and we can use the Thomson scattering cross-section, namely 
A(E e ,T ph ) — > 1. Hence we obtain equation (12a). 
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APPENDIX B 

CONTRIBUTION OF PERTURB ATIVE TERMS IN THE TRANSPORT EQUATION 



Bl. HIGH ENERGY REGION E e > E+ 



Since we can neglect the fluctuation due to the reacceleration in the HE region, we take 
here the second term in equation (25a) alone, omitting the second term in equation (22). 
The transport equation for the electron density in the HE region, Ne°)(r; E e ,t), without the 
perturbative term, is given by 

9 V • D(r; E e )V - e(r)-^W e (E e )} -N^(r; E e , t) = Q(r; E e , t). 



dt 



As discussed in §4.1, we regard Ne°)(r; E e ,t) as the solution of the first order approxi- 
mation for equation (21), so that we have the following equation with the perturbative term, 
n(r)[W n (E e )Ne°) (r; E e , £)]', moving it to the right-hand side, 

jLvJ)(r;£ e )V-e(r)^W e ^ 

Now, we rewrite the solution 

N e , € (r;E e: t) = N<® (r; E e ,t) + N^(r; E e ,t), 



leading to 

d 



d 



dt -V.D(r;E e )V-e(r) W,(E e 



■N^(r;E e ,t)^n(r)-^W n (E e )N^(r-,E e ,t). 



Thus for the steady state (d/dt = 0), we have the solution of the second order approximation 

/•OO 

N^(r;E e )= / J7 e «»(r; y)fi°\y; E e )dy, 
Jo 

with 

'e{r)c^-V-D{r)v\flj?\r-y) = & 0) (r)8(y), 

^^-^>v e (E e )]-/f(y;^) = o, 



where Q (0) (r) and fk } (0;E e ) are given by replacing Q^(r) [= Q(r)] and jf } {0;E e ) (see 
eqs. [24b] and [28b]) with 

g(o)( r ) = g e-^-W/^ =► Q(°)(r) = n(r)e-l»'^^ ^ 

^0 Jk\U k ) 
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/i°)(0; E e ) = E?~* =► /f(0; £ e ) = [Wffifi.) 



with 



J K (u,U)= f fJ K (Ut)dt, 
Jo 



see Table 4 for /c, w K , and £/ K , and J K (U) is the Bessel function of the index with /c (Paper I). 

Corresponding to the replacement of Q^°\r) =>- Q^°\r), the scale heights in the source, 
r Q and z Q , must be replaced as 



1 11 112 

-=*►—, — ^ — + — = — 

Tq Tn %n %D %n 



leading to the following replacements, 



while the radial scale height in the source, r Q , doesn't appear explicitly in this procedure. 
Now the Laplace transform of fn\y; E e ) is immediately given (see eq. [31]) by 

Fr°2(E e ) _ f 00 [W n {Eo)F$(Eo)Y Yr , e (E e ,E ) 



F<$(E e ) J Ee W e (E e )F^(E e ) 

where [•••]' denotes the differential with respect to E , and see § 4.2 for Y r ^(E e) E ), and we 
obtain 

JViS (r; £ e ) 2 J- K (w K , C/ K ) F r ( S (£ e ) 



with 

^ = e r /n r = 77oe- 2r(1/fe - 1/f " } ; r? = e /n . 



(B2) 



In Figure 24a, we show N^(r; E e )/N^(r; E e ) against E e at the SS with [fj Q: f e ] = [2eV, 
8kpc], corresponding to [e Q ,n ] = [2eVcm~ 3 , 1cm -3 ]. Then one finds that the perturbative 
contribution due to the energy change in proportion to the gas density, n(r), is less than 
10% in the energy region E e ^> E+. 

We finally obtain 

N e ,(r; E e ) = ^f^ F^E^e^, (B3) 

with 

F r AE e ) = F^(E e ) + ^f^ F^). (B4) 

'/0 Jk\Uk) 
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B2. LOW ENERGY REGION E e < E+ 

In the LE region, the fluctuation term due to the reacceleration , (AE^ / At) rea , becomes 
now effective as compared to the average energy-loss term due to the synchrotron-IC effect, 
{AE e / At) sic , in proportion to the energy density, e(r). So equation (25b) is approximately 
written as 

n(r)W n (E e ) + 0[e(r)W e (E e )] ~ n(r)[W n (E e ) + fj*W e (E e )], (B5) 

with 

e(r)/n(r) « (e(r)/n(r)) eS = 77*, (B6) 

where the effective value of 77* is of the magnitude of [1-5] eV, and for instance 77* = % ~ 2 eV 
at the SS. 



Neglecting the second term in equation (B5), we have the solution for the principal 
term, corresponding to equation (32), (see §3.3 and Table 4 for u, U u ) 




Electron energy; E e (GeV) Electron energy; E e (GeV) 



Fig. 24. — Numerical results of the contribution of the perturbative terms in the (a) HE 
region and (b) LE region. 
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Now putting 

W:(E e )=W n (E e )+fj*W £ (E e ), 
the electron density with the synctrotron-IC effect (perturbative term here) is immediately 



N*(°)( r -E) - 2 Q° J »( u »> U ^ F* m (E ) e -W ZD 

iv e,n \l i^e) ~ T (TT \ r,n V ' 



(BIO) 



where F*,n\E e ) is given by replacing W n with W* in equations (B8) and (B9), and it is 
related to the perturbvative term, Ne°) , discussed in §4.1 as 



E e ) = JVg(r; E e ) + JV<J>(r; E e ). 



(Bll) 



Here one should be careful of the integral range of E in equation (B8), [.Emm, E max ], since 
we have two zero points in W*(E e ) at two energies, E~ and for instance, [E^,E+] ~ 
[0.1, 1] GeV for £ = 50 mbarn and fj* = 2 eV, and 



[-^minj -^n 



[E e , 


E-] 


for 


Ee < E~, 


[E~ 


E e ] 


for 


E~ <E e < E+, 


[E e , 


+oo] 


for 


E e > E+. 



We present N^} /Nl\n against E e for several sets of [r/*; z n , z D ] at the SS in Figure 24b, 
and one finds it is much less than 10% in the low energy region E e < E+. 

B3. CONTRIBUTION FROM THE FLUCTUATION IN THE REACCELERATION 



Once we confirm that the contribution of 0[e(r)W e (E e )] is approximately given by 
equation (B5) with the energy loss proportional to n(r), it is possible to use the path length 
distribution, IJ n (r;x), as presented in Paper I, 



RJr-x) ~ 2 _^° U *) c -*rx-\z\/z D ^ 

n c J V (U U ) 



but the slab equation is now slightly cumbersome, 



C ^ dE 



—mEe) - <o\^Et a -F r , n (E e ) = cEr 



(B12) 



Remembering that F* J n\E e ) is a solution of the equation 



d 



C*rE: - g^mEe) -E r f\E e ) = cE~\ 
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we can write F rn (E e ) as 

FrAE e ) = F r f\E e ) + F^(E e ), (B13) 



with 



<*rK - ^ W n(E e )} -F^(E e ) = cCo]^E* e -"F r f\E e ). (B14) 



e 



AdE 2 ' 



The solution for equation (B14) is immediately obtained after the following replacement 
in equation (28b) 

1 P2 

and the explicit form is given by 

pW/p \ i r°° r 77i2-a j-i*(0) / r-i 

= CC07 / 1 ° Fr ' n e- y -^^), 

F.t^e) 4 ^ e ° \W*(E e )\F r f\E e ) 

where F r * n (E e , E ) is given by replacing W n (E e ) in equation (B9) with W*(E e ). 

Now, from equations (B12) and (B13), the electron density in the LE region is given by 

N e , n (r; E e ) = ^ ^^ F r , n (E e )e-^, (B15) 

where F rjn (E e ) is given by equation (B13). 

Corresponding to equation (26b), we rewrite equation (B15), dividing into three terms 
as (see also eq. [Bll]) 

N e , n (r; E e ) = N^(r; E e ) + JV£>(r; E e ) + JV«(r; E e ), 

and in Figure 7 in the text, we present N^n/[N^l + N^}\ against E e . One finds that the 
contribution of the fluctuation is effective around GeV region, boosting the electron density 
without the fluctuation by approximately 25%. 



APPENDIX C 

EMISSIVITY of 7's COMING FROM INVERSE COMPTON PROCESS 

In this appendix, we omit the suffix % for simplicity. The production rate of 7's per unit 
time due to the bremsstrahlung, P EB (r; E e , E 7 ) = n(r)c<7 EB (E e , E 7 ) in equation (34), must 
be replaced by 

poo 

P IC (r; E e , E 7 ) = n ph (r; E ph )ca IC (E e , £ 7 ; E ph )dE ph , (CI) 
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where E p h is the energy of the target photon before electron scattering, and E m is given by 
solving equation (A3) with E M = E 7 with respect to E ph (= E m ), 



E m = E m (E e , E y ) = k B T ph 0, 

with 



2 

x i 



n = n (p t \- m ^ 2 / 2 x _ El 



and a IC (E e , Ery] E p h) is the production cross-section of 7's due to IC scattering, which is 
summarized in the right-hand side of Table 5. 

Remarking that the integral range, -E m <-Eph<oo, in equation (CI) corresponds to 
~ (m e c 2 /2E e ) 2 < q < 1, using a parameter q (Blumenthal & Gould 1970) appearing in the 
fourth line of the right-hand side of Table 5, equation (CI) is rewritten as 

P IC (r; E e , Ey) = e ph (r)w T $ IC (x,£ 7 ), 

with e p h in units of eVcm -3 , and see equation (A5) for w T , and 



[Ey/GeV] 



9(1 - -) 2 fw vh {el/ q )<p lc ^,q)qdq, (C2) 

v Jo 



where W p h(k) with k = 0%/q corresponds to the two types of energy spectra of the photon 
gas, the Planck function for the CMB, and the gaussian function for the stellar radiation and 
the re-emission from the dust grains, and 4>i C {x, q) is given in the right-hand side of Table 5. 
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